PHILOSOPHICAL TRANSACTIONS. 



I. Some Measurements of Atmospheric Turbulence. 

By Lewis F. Kichaedson. 
Communicated by Sir Napier Shaw, F.R.S. 

Received October 16, 1919,— Eead February 26, 1920. 

[Plate L] 
Contents. 

Page 
I. Notation 1 

II. Shearing stress deduced from pilot balloon observations . 3 

III. Eddy-viscosity from the same. It varies with height and with direction (anisotropic) . . 5 

IV. Eddy-diffusivity from the dispersal of smoke or of parachutes 5 

Y. General theory of eddy-diffusivity deduced from scattering 8 

VI. Osborne Eeynolbs' eddy-stresses. A measure of gustiness 10 

VII. Summary of theory of scattering of particles of air 15 

VIII. Numerical values from scattering of particles 19 

IX. Cumulus eddies in calm weather 26 

X. Summary 28 

I. Notation. 
The following notation is used throughout. The co-ordinate axes are a right-handed 
rectangular system x 9 y, h, in which Oh is directed vertically upwards, and Ox lies in 
any azimuth which happens to be convenient. Elements of distance to east and to 
north are denoted by de, dn, so that they are special cases of dx, dy. Th6 
atmospheric density is p, the pressure is p, acceleration of gravity is g, latitude <j> is 
reckoned negative in the southern hemisphere, and w is the angular velocity of 
the earth. Velocities are denoted by v with a suffix to indicate the direction towards 
which they blow. Momenta per unit volume are denoted by m x , m^, m H . The 
eddy-diffusivity is denoted by a capital K as in G. I. Taylor's recent papers. 
Another, and in the author s opinion a better, measure of turbulence is £ discussed 
in a previous paper. # The relation K to £ is given by 

dp\ s dp)~~di ~~ ^dh 2 K } 

where x is either potential temperature, or else mass of water or smoke per mass of 
atmosphere. If p and £ were independent of height, then from (l) we should have 

c == Q P -"-• • V / 

* L. F. Kichardson, 'Roy. Soc. Proc.,' A, vol. 96 (1919), pp. 9 to 13. 
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2 ME. LEWIS F. EICHAEDSON ON 

It is suggested that £. might be named " the turbulivity." Its dimensions are : 
(mass) 2 x (length)" 2 x (time)" 5 . 

The advantage of using £ instead of K is that the former enables one to allow for 
variations of density and of turbulence in a simple and natural manner. The 
disadvantage of £ is that it has no name derivable from indoor physics. It is suited 
to the free atmosphere. We might compromise by using in place of K or £ the 
" eddy-conductivity/' c, defined by the equation 

^ = l(4^) or ' app ' roximately 'l = 7^( c l) • • • (3) 

In so doing we gain an acceptable name " conductivity," but we lose by the explicit 
appearance of density in the equation. Either c or £ allows variations of turbulence 
with height to be treated correctly, while K does not do so, as has been pointed out 
elsewhere by the author. # The dimensions of c are (mass) x (length) -1 x (time) -1 . 

This c is of the same dimensions as the measure of turbulence discussed by 
W. Schmidt, of Vienna, under the name of " Austausch " in two important papers. 
('Sitz. Akad. Wiss./ Wien, 1917 and 1918.) 

However much turbulence and density may vary with height 

g 2 p° = £- (4) 

On the contrary if there are no variations with height, 

c = /oK. . o . . . (5) 

The six components of stress are denoted by xx> yy, hh y xy, yh, hx> as in the writings 
of K. Pearson. 

The convention adopted for the signs of eddy-stresses conforms to that of Love's 
" Theory of Elasticity." Tractions are reckoned positive. That is to say, a direct 

stress such as xx is positive if it be a tension, negative if a pressure ; and a shearing 

stress such as xh is positive when the air on that side of a level surface for which h is 
greater (i.e., above), drags the air below in the sense of x increasing. 
The definition of eddy-viscosity adopted in this paper is 

eddy shear ing stress / * \ 

rate of mean shearing strain >•••••••• 

in agreement with the definition used by W. Schmidt (loc. cit. y 1917, p. 5). 

The advantage of this definition is that it is simply based on the fundamental ideas 
of stress and strain, as well as being in harmony with the definition adopted in the 
theory of viscous liquids, (q/1, Lamb, 'Hydrodynamics/ IV. edn., § 326). 

The question may arise as to whether the viscosity defined by (A) can ever become 
infinite by the vanishing of the denominator. The point is discussed by the author 

* Loc. cit. 
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in the paper already cited (p. 13), and the conclusion is reached that such an 
occurrence would be highly improbable. 

The relation of £ to the eddy-viscosity is most easily reached via terms dm x /dt, 
dm Y /dt in the dynamical equations. If pressure-gradient just balanced geostrophic 
wind we should have ^ 

djxh) 3% / 6 \ 

dh dt W 

Now by the definition of viscosity given above 

xh = ju . dv x fdh- (7) 

Substituting (7) in (6) and inserting m x = pv x there results 

idhV-dh\-~dT (8) 

It is seen that this equation becomes identical with (l) if 

x = v* and i = g 2 P v, . ■ (9) 

of which the latter is the required relation. 

On comparing equations (8) and (3), it is seen that eddy-viscosity, /x, and eddy- 
conductivity, c, are of the same dimensions, and appear in their respective differential 
equations in the same way. Indeed, Taylor has suggested that they are equal. # 
This likeness would be a good argument for recording observations in terms of these 
two quantities instead of in terms of diffusivity K or turbulivity £ 

II. Shearing Stress from Pilot Balloon Observations. 

(Condensed and revised January 22, 1920.) EKMANt in a remarkable paper 
pointed out that the total momentum of water produced by a tangential stress on 
the surface of the sea, in the steady state, is directed at right angles to the tangential 
stress, and its amount is quite independent of the value of the viscosity or of the 
variation of viscosity with depth. The same applies to the atmosphere. We may 
use this principle to find the shearing stress on the ground, provided we have a 
measure of what the momentum would be if the surface stress were zero. 

I have taken the wind at a height of lj km. to 2J- km. as the standard of 
reference, because, by so doing, the term depending on curvature of path, and other 
small terms in the dynamical equations, are automatically allowed for to a first 
approximation. The stress at 2 km. is undoubtedly much less than that on the 
ground, and is neglected. It is best to select observations in which the momentum 
becomes nearly independent of height above 1*5 km. A table of results follows. 
They were computed with the help of Mrs. L. F. Eiohardson. Dr. H. Jefferys 
says the selection will select abnormal lapse-rates and so abnormal viscosities. 

* < Phil. Trans./ A, vol. 215, p. 22. 

t " On the Influence of the Earth's Eotation on Ocean Currents," by V. W. Ekman, s Arkiv for Matem. 
Astr. och Fysik,' Stockholm, Bd. II., No. 11 (1905). 
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SOME MEASUREMENTS OF ATMOSPHERIC TURBULENCE. 5 

III. Eddy- viscosity from Pilot Balloon Observations. 

(Abridged, January 22, 1920) The method of the last section will give the 
difference of the shearing stresses on two surfaces of any horizontal slab of air. If we 
choose one of the surfaces so that dvjdh = and consequently also the stress 
vanishes, we obtain the stress on the other surface. This has been done for some 
very smooth means for Lindenberg (E. Gold, Met. Office, ' Geophys. Mem./ V., p. 143). 
The results are set out in the following table : — 

x .ajrle 1 jl. ■■" '""xjinden oer g . 

The x axis is directed with the surface wind. The stress is that exerted by the 
upper on the lower layer. 



Height above 


Eddy-shearing-stresses. 


Rates of mean shearing. 


Eddy viscosities. 






«? 




Parallel 


Perpendicular 


mean sea, 






flV*. 


riiL, 


to wind. 


to wind. 


kilometres. 


xh 


yh 


dh 


KJUy 

dh 


xh 


yh 




dynes cm." 2 . 


dynes cm."" 2 . 


sec." 1 10 3 x 


sec." 1 10 3 x 


dv x Jdh 
dyne cm." 2 sec. 


dvy/dh 
dyne cm." 2 sec. 


1-0 


+ 0*05 


-0*09 


- 1*0 


- 0*2 


50 


450 


0-8 


+ 0*01 


-0*37 


- 1*2 


- 1*5 


10 


250 


0-7 


zero 












0-6 


+ 0*02 


-0*71 


1 


- 3*2 


20 


' 220 


0*4 


+ 0-09 


-0*99 


6*0 


- 9*0 


15 


110 


0'3 


+ 0*39 


-1*06 


12*5 


-12*2 


31 


87 


0-2 


4-0-69 


-0-90 










0*12 ground 


+ 1-10 


-0*64 


21*51 


-17*51 


51 


37 



To obtain a quantity comparable with £ we must multiply the eddy-viscosity by 
g 2 p which is approximately 1100 c.g.s. units. 

The mean of the viscosities in the two directions increases with height as we 
might expect from other observations (vide Part VIII. , below, also 'Roy. Soc. Proc./ 
A, vol. 96 (1919), p. 18). But the most interesting thing about this table is the 
marked lack of isotropy in viscosity. The air appears to he more viscous, for large 
motions, across the wind than parallel to it, except just near the ground. 



IV. Eddy-diffusivity from Smoke or Floating Bodies. 

Some direct measurements have been made by observing the gradually increasing 
scatter of smoke or other visible material carried along by the air. The changes in 
height of a large number of small portions of air are observed during a fixed interval 
of time. These changes are found to be distributed about their mean value 
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approximately according to the ordinary "law of error." Their scatter in height is 
measured by the "standard deviation/' computed by the familiar methods. # In 
order to find the diffusivity K the observations are compared with an appropriate 
integral of the approximate equation 

*=K^ (1) 



dt dk 



A i e nu+A,) (2) 

a/4£ + A 2 



Such a one is 

X = 

where A 1? A 2 , A 3 are constants. 

This integral represents a horizontal lamina in which the density x is distributed 
about a mean height A 3 according to the law of error. The square of the standard 
deviation of the mass in the lamina can be shown to be 

(4* + A 2 )| (3) 

So if the scatter of the same set of particles be observed at the beginning and at the 
end of an interval T of time, it follows that 

K = -= (increase during T of square of standard deviation). . . (4) 

But the increase of the square of the standard deviation is equal to the square of the 
standard deviation of the change of height. Accordingly 

K == — (square of standard deviation of change of height during T). . (5) 

Jj JL 

In this last transformation we have assumed that K is sensibly independent of 
height. This is permissible because the range of scatter can usually be made small. 
For the same reason the density of the air may be taken as independent of the 
height, so that we may obtain from K, the constant £. which we require when 
pressure is taken as independent variable in place of height, in accordance with (l) 
above. This procedure is not perfectly satisfactory but it is very convenient. It 
gives £ = # 2 /> 2 jK and " eddy-conductivity " = />K. 

There is no need for the changes in height to be simultaneous for all the portions 
of air, and in practice it is much more convenient to let them be successive. 

Varieties of particles. — I have observed the scattering of smoke from a 
smouldering wick, from burning weeds, from factory chimneys and from ship's 
funnels : also the scattering of portions of cloud near the horizon and of puffs of 
ammonium chloride from a special apparatus. Chimney smoke is not to be 
recommended, as it rises through the air. Clouds and steam may mislead one by 

* Vide < Computer's Handbook,' M.O. 223, Section V. 
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evaporating. The cold NH 4 CL smoke proved more satisfactory in these ways. 
Lycopodium dust might be better still, as isolated grains would fall at a definite 
rate relative to the air. But I have not succeeded in making the lumps break up 
into grains. The downy parachute which carries the seed of the dandelion, Taraxacum 
officinale, has been found to be convenient. When the seed is broken off, the 
parachute falls at a rate of 10 to 15 cm. sec." 1 in still air. The standard deviation of 
this rate must be allowed for. The formulae for correction are given below. A brown 
parachute, twice as large each way as that of taraxacum, grows near Benson. I am 
indebted to the Botanical Department of the British Museum for a search among 
their dried specimens for a large white parachute. A splendid one came from an 
African plant called strophanthus. The parachute is about 6 cm. in diameter and has 
a long stalk by which it can be held conveniently. When the seeds were broken off 
the parachutes fell, in still air, at an average rate of 20 cm. sec." 1 . 

Other artificial clouds, which have been used with success, are paraffin-oil vapour 
from an extinguished blast lamp, and smoke of phosphorus pentoxide made by 
dropping calcium phosphide into dilute hydrochloric acid. In strong winds the smoke 
from a firework known as " Vesuvius " # is convenient. 

If one could mark and follow individual molecules equation (5) would give the 
molecular diffusivity in still air, 0'2 cm. 2 sec." 1 . Actually what we observe is the 
centre of a small puff of smoke, and this is not constantly the position of the same 
molecules, so that in still air we find K = 0. To be perfectly exact all observations 
of K by this method should be increased by 0'2 cm. 2 sec." 1 , an entirely negligible 
correction. In any theory the diffusivity depends on the motions which the theory 
does not follow in detail. In laboratory experiments, in which the molecular motion 
only is ignored, K is taken as 0'2 cm. 2 sec." 1 . In meteorological telegraphy variations 
of wind of less than 10 minutes' duration are ordinarily ignored, and there is an 
appropriate, much larger, value of the diffusivity. In a certain scheme for numerical 
prediction it is proposed to average the wind over periods of 6 hours, and the 
further variations thus omitted must be taken into account by further increase in K. 
It follows that the puffs of smoke should be so small as to allow the smallest eddies 
to be observed, and, for the last-named purpose, that the observations should be 
spread over a period of 6 hours. In obtaining the data in the following table 
I believe the former condition has been fulfilled, but the latter has not. When only 
the order of K or £ is required, it is enough to assume that the standard deviation 
is | of the distance between the extremes of height observed, when the number of 
observations is about 40. 

Observations very near the earth's surface have peculiarities. It is obvious that 
dxfdh = at an impermeable horizontal surface. This condition can be satisfied in 
the integral by taking the portion of the distribution which would be cut off by the 
surface, reflecting it in the surface, and adding it to the rest of the distribution. 

* Made by Messrs, C. T. Brock & Co, 



8 ME. LEWIS F. KIOHAEDSON ON 

The standard deviation is not then sufficient to give K. It is necessary to make 
more elaborate computations with K. Pearson's "incomplete normal moment 
functions." 

These difficulties are avoided in the symmetrical case when the source of smoke is 
exactly on the ground, and the smoke does not rise or fall by its temperature. 

The results of these measurements are set out in the last column of Table IV. It 
is seen there that K increases from 5 near the surface of land up to 10,000 at the 
height of a factory chimney. 

But before considering the results further a fuller mathematical investigation will 
now be made. 

V. General Theory of Eddy- diffusivity Deduced from Scattering. 

The foregoing theory of the diffusion of a lamina assumes that the diffusivity is 
constant throughout the space, and that the density in the lamina does not vary 
except in the smooth regular manner indicated by the " law of error/' But it is well 
known that the wind has an intricate structure. Thus if observations of the smoke 
puffs are to yield a measure of the diffusivity from the formula 

tj* • -± increase in square of standard deviation 

dinusivity = - — -. - 1 — ^ : -. — ~-. — , 

twice corresponding increase in time 

then either the interval of time in the denominator must be long compared with the 
fluctuations of the wind in time, or the initial standard deviation must be large 
compared with the fluctuations of the wind in space, or both conditions must hold. 
The former condition is an inconvenient one in practice, because puffs are apt to fade 
before a sufficient time has passed. Dandelion parachutes, with the seeds removed, 
may be better than smoke for this purpose. 

The following theory brings to light some of the assumptions involved in the 
measurement of diffusivity by smoke puffs. It was contrived . specially in order 
to avoid "the distance through which an eddy moves before mixing with its 
surroundings," a quantity which occurs in Taylor's theory, but which does not lend 
itself easily to measurement, except in the case of cumulus eddies. See Section IX. 
below. 

The potential* temperature 6 does not change at a point moving with fluid, if 
radiation and precipitation can be neglected. Now let a portion of an eddy move from 
a height h x at time t x to a height h 2 at t 2 . Then, regarding 9 as a function of h and t, 
we have 

o {K, h) = e {h* t 9 ). (i) 

* Potential temperature is the temperature which the air would acquire if compressed adiabatieally 
to a standard pressure. If 6 is to be of service in dealing with cloudy air the standard pressure must be 
high enough to evaporate the cloud in all samples. 
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From each side of (l) subtract G (h l9 t 2 ) and divide through by t 2 —t x . Then 

e (h l2 Q-e (K t 2 ) = e (h 2 , Q -e ( K t 2 ) (2) 

i 2 ^ 1 i 2 v ± 

Now the left side of (2) is the finite difference ratio y at the height h l9 which is 

what we want to find in terms of the spacial distribution of 9. Expand the right- 
hand side of (2) in powers of h 2 — h x by the well-known theorem in the calculus. It 
follows that, if subscripts indicate the time and height 

^\^ = ^{®j^- K H^L iI ^ + ^ hwi ^ ■ ■ (3) 

The difference ratio on the left of this equation is centred at the same height h x as 

the differential coefficients on the right of the same, but at a time jf(t 2 —ti) previous. 

This slight misfit in centering will not matter, because t 2 —t 1 will be of the order of 

» 

one minute or less, whereas we are next going to take the average of each term in 
(3) over a much longer time, say 6 hours. The subscripts may now be omitted as 
unnecessary. Let a bar over a symbol, or group of symbols, denote the mean value 
over this longer period. Let a dash denote the instantaneous deviation from this 
mean, so that, for instance, we have for every fluctuating quantity a formula such as 



dh \dh) + \zh) ' ' m : • ( 4 ) 

Now the mean of any dashed quantity vanishes. (5) 

Again the mean of the product of any dashed quantity into any barred quantity 

also vanishes.- (6) 

We shall further suppose that the mean of h 2 —h x vanishes, . . .. . . . (7) 

that is to say that there is no mean vertical displacement. 
Then, in the first term on the right of (3) 



39 
dh 



(th-K) ={(1)+(1)} {(^-M+fo-Aiy} 



^ (A 2 — ^i) (8) 



dh 
because of (6) and (7). 

/r\f\\ f 

Now i — )(A a — Ai) ; , after being divided by t 2 —t x and by the " standard deviations" 



of(— ) and of (h 2 — h x )\ becomes equal to the correlation between dQ/dh and 
vol. ccxxi.— a. o 
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{h 2 —h-)l(t 2 —ti). So the first order term on the right of (3) vanishes, on taking the 
mean, if the variations of oQ/dh, in time at a fixed point, are not correlated with the 
variations in (h 2 —h 1 )/(t 2 —t 1 ), at the same point and time. 

In cumulus cloud eddies, the variations of velocity are caused by variations of the 
potential temperature 9, so that a correlation is almost certain to exist. On the 
contrary, when the eddies are due to dynamical instability, the correlation may be 
expected to vanish. In the latter case, it is the second order term of the right of (3) 
which becomes effective, so that _ 

W-WWzEl (a) 

*8t ~dh 2 '2(t 2 -t 1 ) ••••••••• W 

Now suppose further that either 3 2 6/3A s has no variations at a fixed time and level, 
or else that its variations are not correlated with those of (h 2 —h 1 ) 2 . Then (9) 
simplifies to 



SO 3 



£%=*# (10) 



%t dh 2 ' 2(« 8 -0 
Thus 

.......... _ * 

(h 2 — h^fj2 (^ 2 ~0 is the eddy-diffusivity K (ll) 

It is seen to be identical with that derived in Part IV. above, by considering the 
diffusion of a lamina, in which the density was distributed according to the law of 
error. It is a quantity easily measured. 

Of course if t 2 -~t 1 were sufficiently small, say T V second, then it would be the first 
power of h 2 — h u which would be proportional to t 2 —t 1 , instead of the square. This 
suggests that t 2 — t x must be long compared with the fluctuations of the wind. On 
the other hand t 2 —l 1 must be short compared with the period, of say 6 hours, over 
which the averages denoted by the bar are desired to be taken. 

A similar argument can be applied to any other quantity which, like 0, does not 
change following the motion of the fluid, provided it has space-rates independent of 
the time-variations of velocity. Thus the mass-of-water-per-unit-mass-of-atmosphere 
may replace in (10) with similar restrictions. 

When we consider diffusion in three dimensions there may be six coefficients of 
diffusivity corresponding to the six components of stress. 

VT. — Osborne Reynolds' Eddy-stresses. 

But we cannot, without further investigation, apply the preceding argument to 
the diffusion of horizontal velocity in a fixed azimuth. 

Something might perhaps be deduced from the well-known theorem that, when p is 
constant, 

where yjr is the gravity potential. 
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The eddy-viscosity can however be measured rigorously ,J3y smoke-puff observations 
made in such a manner as to fit in with Osborne Reynolds' theory of eddy-stresses. # 
This theory is remarkably free from assumptions which might limit its generality. 
It is to be found in Lamb's ' Hydrodynamics/ IV. edn., Art. 369. 

The equations of motion are three, such as 

+ P —- — 2oo sin <pv Y — c (i—^ HVVjj .... (13) 

where c is the " molecular" or u ordinary" viscosity. Note that there is no need to 
assume p to be independent of position. Reynolds assumed this, but for a reason that 
does not concern us. It will be necessary however to assume that //, the variation of 
density at a fixed point, is so much smaller in comparison with p than is v r in 
comparison with v 9 that we may put p = 0. This being so, we find on taking the 
mean that (.13) becomes 



^} - !+'£-*- si " *- ° (* -*" + N 



+ g" (pvx • % + Xx • v' x )+ — (pv x . v Y + pv' x . v' Y )' 



+ -r (p%.% + M) • (14) 

d/fc 

The left side of (14) is the difference between pv x at the beginning and at the end of 
the period through which the average is taken, divided by the period ; and that is 
what we want. The right side of (14) is of exactly the same form in the mean 
quantities jo, v X9 v Y , v R as (13) was in the corresponding instantaneous quantities 
P> v x> v y> V n I except that there is added a force per unit volume in the x direction 
equal to minus 

<y *\ ( p\ ______ 

g-(/^Wl)+ g- (/^Vl) + g^ UVh). ...... (15) 

On working out the corresponding equations for the y and h components, it is seen 
that this additional force per unit volume is just that which would be given by the 
following systems of stresses 



~> 



xx = -pv' x v' x ; yy = -pv f Y v f Y ; hh = -^V H 



> . . . (16) 



xy = -pv' x v' Y ; yh = ~/^Wh ; hx = -^Wx 



* Major G. I. Taylor tells me that he attempted to measure xh with a balloon on an elastic tether 
in 1914. 

c 2 



PlHE 
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when any symbol such as xy is the force in the x-direction per unit area of a plane 
normal to the y axis. Tractions are reckoned positive, as usual. The above is taken 
from Osborne Keynolds' theory, adapted and slightly generalized to suit our needs 
for a rotating atmosphere, having density diminishing with height and a molecular 
viscosity which is not neglected. 

One may form a clear mental picture of these eddy-stresses by imagining the 
scattering of smoke puffs. Let a puff emerge from a pipe at the origin of the 

co-ordinates. After a short interval of time r, let the 

puff appear at the point P on the diagram, as seen by an 

observer at a distant point on the ?/-axis. Now let the 

observation be repeated for a large number of puffs in 

v * • succession, the time r being kept the same for each. 

.; ' ; \ ( We thus obtain a diagram, with a large number of points 

. ; ■ on it, showing the scattering of the puffs after r. Then 

the eddy-stresses are simply related to the correlations 

and standard deviations of this scatter-diagram — under 

certain conditions. For let X, Y, H now mean the 

co-ordinates of any one of the dots on the diagram reckoned from the source of 

smoke. 

Then the velocities of the corresponding puff were 

"Y" TT 

v* = — ; '"h = — , (17) 

provided the time r was so short that, during it, the velocity may be regarded as 
uniform and in a straight line (18) 

Again, the velocity of the puff will be equal to that of the air which it has replaced 
provided the puff is at the same temperature as the air, and provided that the pipe 
points parallel to the Y-axis so that the impulse with which the puff leaves the pipe 
does not show in the projection on the plane XOH. 

Let us suppose that a number of puffs, n in all, are observed. In order to 
correspond with the time-mean taken over 6 hours, which was used in deriving the 
eddy-stresses from the equations of motion, these n puffs should be spread uniformly 
over a similar interval. 

From the scatter diagram we can compute first the mean velocities. For the 
mean velocities are 

v x = 1 = -i- SX ; v Y = T/r = — , • • • • (19) 

r r . n r . n 

where S has the meaning : — -take the sum of what follows it, for n puffs. 
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Then the normal eddy-stresses, xx and hh, are found thus 



xx = - P v' x v' x = -p{v x -v x ){v x -v x ) 



P 



r'n 
So 



2(X-X)(X-X) (20) 



tKJ*AJ 



= —pvx 2 fa\ (21) 



where cr x is the "standard deviation" of the dots in the ^-direction. If o- H is the 
corresponding quantity vertically 

So the direct eddy-stress in the direction of the wind is intimately related to the 
gustiness shown by a tube-anemometer. 
The shearing eddy-stress 

xh = ~"p^ x V H =-4^2 (X-X) (H-H). (23) 

So 

S=-pT-V xn< r x , n , (24) 

where r XH is the correlation between the co-ordinates X and H of the dots. 

By projecting the puffs on the other two co-ordinate planes we should be able to 
measure similarly the remaining components of eddy-stress. 

To find the eddy- viscosity we must compare the sheering eddy-stress xh with 
1^ + ~3/ ? w ^ c ^ is. the rate of shearing strain in the mean motion. Usually dv B /dx 

Ofb OX J 

is negligible, so that the rate of shearing is dv x fdh, a quantity which can easily 
be observed. At first sight one might think that dv x /dh was simply related to the 
slope of the regression line in the scatter diagram ; but on examination this proves 
not to be the case. The slope of the regression line is independent of t, because (18) 
is satisfied for all permissible intervals of time. 

It should be noted that no shearing stress such as —pv' x v'n can exceed, in absolute 
value, the geometric mean of the corresponding pair of direct stresses —pv f x v r x ^ 
—pv'nv'n for the same reason that a correlation coefficient cannot exceed unity. 

The probable errors of eddy- stresses, determined from the scattering of particles 
moving with the air, may be taken to be as follows 

Probable error of xx = 0*674 . xx a/ — ........ (25) 



<*" "■"» / /y» fa _j_ /y*/y* fa fa 

Probable error of xh = 0'674 ^ / tt— * — - (26) 
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These are based on the assumption that the number of particles is not less than 
say 20, and that the scatter is "normal," so that the logarithm of the density of 
particles would be a quadratic function of x, y, h if a very large number of particles 
were observed. Reference may be made to the fundamental papers on probable 
errors: Fiixm and Pearson, 'Phil. Trans./ A, vol. 191; W. F. Sheppakd 
' Phil Trans./ A, vol. 192. Equation (26) follows from the probable error of the quantity 
called by K. Pearson the "product moment coefficient taken about the mean" 
Equation (25) follows from (26) on putting x == h, or may be deduced independently 
from the probable error of a standard deviation taken about the mean. 

Corrections for the Motion of the Parachute Relative to the Air. 

When observing eddy-stresses by the aid of the parachutes of plant seeds it is 

desirable to allow for the velocity of the parachute in still air. For large specimens 

of the parachute of Taraxacum officinale, after cutting off the seed, the velocity in 

still air was found to have a mean of 12 cm. sec." 1 with a standard deviation of 

2 cm. sec." 1 . It may be shown from the equation of motion that this limiting velocity 

is acquired in a negligibly short interval of time. Thus call the limiting velocity 

c downwards, the instantaneous velocity downwards u. Then if the friction is 

proportional to the velocity 

g (mass) = c x F, (27) 

where F is a constant. But, when accelerating, 

(mass) -. — — <7(mass) — uF: (28) 

Eliminate the mass between these two equations and there results 

| *<«=£> = - <„-,). . ........ (29 ) 

So that the discrepancy between the actual velocity u and the terminal velocity c 
sinks to e" 1 of itself in a time equal to c/g, which for the taraxacum parachute having 
c = 12 cm. sec." 1 would be only a hundredth of a second. 

Less negligible is the variation of the velocity c from one parachute to another. 
What we actually observe is not v K the upward velocity of the air, but v R — c. Now 
write V for the mean velocity of the parachute in still air, and d for the deviation from 

the mean. Then in finding the direct stress hh we must perforce work out first the 

"raw" moment {(v K — c) 2 } = [{(^ H — c) + (d K — d)Y\ On expanding and remembering 
that a bar put over the product of a dashed and a barred symbol, causes the result to 

vanish, and also that v' a —d = 0, it is found that 
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it follows that the corrected stress is given by 



to =~ P iV H =- P {( % - c )(%~c)-(i" n ~c) 2 -cV} (30) 



which is the actual formula used in working up the observations. Here cV is the 
square of the standard deviation of the velocity of the parachutes in still air. 

Next we require the product moments in order to find the shearing stresses xh and 

yh. The raw moment v x (v n —c) is found on expanding to be equal to v x {v K —c)-{-v f x v f 1I . 
So that the corrected value for the stress is 



xh = -pVx = ~p {v x (v K -c)-V x (v K -c)} 9 ..... (31) 



and as cV does not appear, the scatter of the velocities of the parachutes in still air 
does not make a correction necessary for the shearing stresses. 

VII. Summary of Theory of Scattering of Particles of Air. 

The conclusion we have reached is the following. For any sort of eddy, whether 
due to " dynamical instability," or to the rising of heated air in cumuli, the eddy- 
stresses are best measured by equations (22), (24) and the like, because the theory 
from which they are derived is very general ; and the eddy- viscosity is best measured 
as the ratio of the shearing eddy-stress to the rate of mean shearing strain. It is 

conceivable that xh found from (24) might turn out to be zero. In that case it would 
be necessary to investigate effects of higher order. This might possibly be done by 
developing, for the quantity (^pif + py/s+p) in equation (12) an analysis similar to 
that of (l) to (11) for potential temperature. The diffusivity for potential 
temperature, on the other hand, should be measured differently according as the 
eddies are produced by variations of potential temperature or not. Thus for cumulus 
eddies we should take the mean of (3), retain the linear term on its right-hand side 
and neglect the quadratic one. Then t 2 —t x must be small, so that 



m _ 36 (h 2 -Jh) „ /30V / 



«* ' a/i* t 2 -t x \dhj 



( -r~r ) v\i> (3lA) 



The diffusivity is measured as the right side of this equation divided by 3 2 9/3/?, 2 . 

nnr 

X I1US /^/^\\ t 

k _U Vh (32) 

d 2 Q/dh 2 

But for eddies due to dynamical instability, neglect the linear term in (3) and measure 

the diffusivity as m ■» \2 

\'h 'Hi /^^ 

I V, • . . . . • • • • • \OtJt 

2(« a -0 

where (£ 2 — O must not be too small. It will be interesting to see whether eddy- 
diffusivity is found to be equal to eddy-viscosity divided by density. 
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Let us now see how these two modes of disposal of smoke fit together — -the one. at 
short intervals of time, governed by the eddy-stresses, the other at long intervals 
governed by the diffusivitjr. Let us draw the trail of smoke as it would appear to a 
distant observer looking, say, horizontally. The scatter may be represented by 
drawing, two lines through the points where the smoke-density has its standard 
deviation in height above or below its mean. Let T be the time since the smoke 

emerged from its source. Suppose that the stress hh 9 the diffusivity K and the mean 
velocity v x are all constant along the path of the smoke. The time T taken to travel 
a horizontal distance x measured down the trail from the source is x/v x . Near the 
origin T = t in equation (22), and so the standard deviation in height is 

*i<^w, (34) 

representing a pair of straight lines intersecting at the source. Further down the 
trail T = U—t Y in equation (ll), and the standard deviation in height is 



± vx 



. /\ / — • • • « > • • • • • \ / 

y v x 



representing a parabola with, horizontal axis and its apex at the source. Thus 
according to the theory, the smoke may be said to be contained within a paraboloid 
which has had its blunt end sharpened into a cone. The preceding theory gives us 
no clue as to the manner of transition from the cone to the paraboloid. 

When we can observe a sufficient length of the path traced out in space by a single 
small portion of air the eddy-stresses and the eddy-diffusivity may be deduced from 
the irregularities in the motion. With this object I have observed the motion of anti- 
aircraft shell bursts, and of portions of cloud, by means of an Abney level or a pocket 
sextant. With better instruments this method might yield a good deal of information 
about eddies at heights such as 2 to 5 km. One principal difficulty is that the shell- 
burst fades away after about 5 minutes, before a sufficient length of path has been 
observed to give the diffusivity. 

If the path is sufficiently high and long, the hills, trees and houses on the earth 
may be regarded as blending into a " roughness." Suppose this roughness to be 
uniform. Then if we had been causing smoke to issue in puffs from a fixed pipe, we 
should presumably have obtained the same scatter diagram for the puffs, within the 
limits of probable error, at whatever point of the path we had placed the pipe, or at 
whatever time we had begun to observe, within limits. If this is so, we may form 
the scatter diagram by taking its origin at every point in succession of the trajectory 
of the single particle. For instance, Captain Cave in his book on "The Structure of 
the Atmosphere," gives several diagrams of the irregularities of height of a balloon 
observed by two theodolites, the uniform vertical motion of the balloon relative to the 
air having been eliminated. From his figure 30, of an ascent on February 19, 1909, 
the following has been deduced, by taking the origin of scattering at every available 
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minute mark, and working out the standard deviation, in height about the mean 
height, after a time a. 

Table III. 



a, sees. 


1x60 


3x60 
0-82 


5x60 


7x60 
0-55 


9x60 


cr 2 

^ - 10 4 X 

2a 
cm. 2 sec." 1 


I 0-93 


0-55 


0-33 


OVr 2 

dyne cm," 2 


I 0-345 


o-ioo 


0-041 


0-029 


0-014 


Number of points 


14 


12 


10 


8 


6 



Now <r H 2 /2a would be the diffusivity if a were a " long " time ; and the criterion of 
sufficiency in length is that <r R 2 f2 a should not vary with a. The small number of 
points makes the probable errors large, so that the decrease of cr H 2 /2a between 
a = 7 mins. and a = 9 mins. is not significant. It looks as though the diffusivity K 
were here of the order of 0'5 x 10 4 cm. 2 sec. -1 . 

Again — pcr H 2 /a 2 would be the eddy-stress, hh, if a were so small that further decrease 
made no further change in the quantity. This stage is not reached at a = 1 minute. 
All that we can say is that the eddy-stress is probably greater than 0*345 dynes cm." 2 . 

The mean height of this observation is 1 km. above ground, and the mean velocity 
1270 cm. sec."" 1 . 

The photograph of a smoke trail in fig. 2 suggests that the eddying is partly 
random but also partly sinusoidal. Let us therefore see what would happen if the 
path of the particle were an exact sine curve without any random variations. Let 
the height h of the particle be given by h— B = A sin qt where A, B and q are 
constants. The increase in height in *a time a would be 



A 



sm^^+|j|--sin|gU--? 






= A . 2 cos qt . sin ( q 



a 



by trigonometry. So that the standard deviation <r H after a time a would be given by 



4A2 sin f / r-* 

^ = \ Al (cos qtfdt 

L J t - o 



where L is a very long time. The integral is equal to |L plus a negligible oscillatory 



2 



It follows that the stress hh, which is the 



part. Consequently o- H 2 = 2 A 2 ( sin ^ 

limit of —p(T K 2 /a 2 when a is small, comes to — pq 2 A 2 /2 ; whereas the diffusivity K ? 
which is the limit of <r H 2 /2a when a is long, comes to gero, 
YOL. coxxi.— a, P 
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Comparison with Taylor's expression for the Diffusivity and with W. Schmidt's 

" Austausch" 

In Taylor's remarkable investigation (' Phil. Trans./ A, vol. 215) from which the 
present research took its stimulus, the diffusivity K is given, in the present notation, 
as the mean value of v R (h—h ) over a large horizontal plane; and it is stated that 
h—h is the height through which an eddy moves from the layer at which it was at 
the same temperature as its surroundings, to the layer with which it mixes. This 
definition of h — h is puzzling, for it seems impossible to reconcile the supposed 
starting and stopping of the air, with the ceaseless motion which we observe in nature, 
except in the case of cumulus eddies. Happily we are now in a position to clear away 
the mystery. For it has been shown independently in the present paper that the 
diffusivity is given by 

K = (h 2 -h 1 y/2(t 2 -t 1 ) 

where (£ 2 —"0 m a time long compared to the fluctuations in the wind, and where the 
bar implies an average taken over a still much longer time. As {t 2 —t^j is the same 
for all the quantities which are averaged, we may remove it from under the bar, 
writing 

K(* 3 -^=i(V^. 

Differentiate this equation with respect to t Jy 



$ 

K = (h 2 -h ± ) -=- {h 2 -h x ) = (A a -^i) • ^h at t 2> 



dt 

thus K is expressed as the mean of the product' of the rise in height during a long 
time into the vertical velocity at the end of that time. It may also be taken at the 
beginning. Comparing with Taylor's form quoted above we see a strong resemblance, 
and we are led to suppose that Taylor's theory makes two unnecessary and 
unnatural restrictions : (l) that the portion of air should start at the same temperature 
as its surroundings ; (2) that the portion of air should finally mix with its surroundings. 
But that if these restrictions be removed, then another becomes necessary, namely 
that (t 2 — tj should be sufficiently long (several minutes). Whether the average be 
taken over a large horizontal plane, or over a very long time (6 hours), appears to be 
a matter of indifference. 

The extent to which Taylor assumes viscosity to be independent of height in his 
general theory ('Phil. Trans./ A, vol. 215, pp. 11 to 13) is this: he neglects the terms 
due to the initial eddying in his equation (6). That is a doubtful proceeding, unless 
the initial eddying is zero : but zero is independent of height. 

The " Austausch " of W. Schmidt is defined by him (in 'Sitz. Akad. Wiss./ Wien 
(1917), pp. 4 to 5) as 

2 (element of mass crossing horizontal plane) x (vertical displacement of element) 

(whole area) x (time of motion) 
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If we replace the summation by an integration over a large area S in the plane of xy, 

the element of mass crossing per unit time is dx dy pv R , so that the " Austausch " 

1 ff ' 

becomes - \pVn (h—h ) dx dy which is p times Taylor's diffusivity as defined in ' Phil. 

b J J 

Trans./ A, vol. 215, p. 3. Thus we must suppose that W. Schmidt's definition of 
" Austausch " requires amplification concerning the interval of time and concerning the 
position in it of the velocity, just as Taylor's definition of K does. 

VIII. Numerical Values Derived from the Scattering of Particles. 

Fig. 1 is a photograph^ of the trail of paraffin vapour from an extinguished blast- 
lamp which projected the vapour in a direction at right angles to the wind. It 
shows a cone, with a blunt point due to the finite size of the source of smoke, passing 
smoothly into a form, which certainly diverges less rapidly than the initial cone, and 
which looks like a paraboloid. Opinion might differ slightly as to where to draw the 
lines corresponding to the standard deviation of smoke. In a " normal " distribution 
0*68 of the whole number of particles lie between the two standard deviations. If 
the lines are placed as in the accompanying black and white drawing, then it follows, 
as the mean velocity of the smoke was 17 metres/sec, and the density of the air 
was 1*21 x 10~ 3 grm. cm," 3 , that 

stress hh = —0*73 dyne cm."" 2 diffusivity = K = 240 cm. 2 sec." 1 

turbulivity = £ = 340 grm. 2 cm."" 2 sec.~ a . , 

This photograph was taken in the evening, when the day-wind was diminishing. 
The source was 190 cm. above ground. Obstructions to windward only subtended 
an angle of 2°*1 at the source of smoke. The exposure lasted 60 seconds. 

Fig. 2 was taken five minutes later in the same place, with an exposure of 
85 seconds. 

The velocity of the smoke had decreased to 1 *3 metres per sec. The measurements 
yield 

hh = —1*2 dynes cm." 3 ; K = 750 cm. 2 sec. -1 ; £ = 1050 grm. 2 cm."" 2 sec." 5 . 

In this case the photograph shows a distinct neck between the cone and the paraboloid, 
at a distance from the source roughly 1'3 times its height above ground. This neck 
can also be recognized in some other photographs. Its presence signifies that 
the motion of the air was compounded of (i.) a random eddying, plus (ii.) a wave 
motion in which the particle of air executed a wave having a length, relative 
to a point fixed to the earth, roughly 2 '6 times the height of the particle above 
ground. 



* Taken at Benson. 
D 2 
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Table IV. — Observed Values of 

The co-ordinate axes are taken to point : Ox horizontally with the mean wind at 
quantities are in G.G.8. units. Tractions are reckoned positive, so that xh and yh 



Entropy gradient 
vertically upwards. 



Surface. 



Methods. 



Positive sunrise . 



Tall grass . . . 



Smoking wick 



Small ? 



Oat field . 



NH 4 C1 puffs. . 
NH4CI puffs . 
Dandelion down 



Small ? 



Small 1 



Corn field . 



Moor . . . . 



Moor, trees 



Smoke puffs . 



« • 



Smoking wick . . . 



• . 1 



Thistle down 



. . 



Positive ? 



Flat field . 



• • 



1919, Mar. 4 



1919, Apr. 7, 
18h. 



o 

r— < 

a> 

r-fl 
4-3 



Sea 



- ^ i 



1919, July 4, f gg 
llh. over 
cast 






Paraffin vapour . . 



Steamer's smoke . 



Sea 



1919, July 5, 
12h. 






Sea 



Steamer's smoke . . 



Steamer's smoke . . 



Sea 



i 



Steamer's smoke . . 



Ditcham . . . 



Flat fields . . 



Capt. Cave's balloon . 



Phosphorus pentoxide 



Notes. 



h 

cm. 

above 

surface. 



I. 



II. 



• • 



III. 



IV. 



V. 



VI. 

Figs. 1 
and 2 



p. 16 



f VI 
1 Kg: 



VI. 

; 3 



cm. 
sec. 



30 

170 

300,000 



10 
60 

200 



dynes per 
square 
centi- 
metre. 



40 

100 

50 



120 



20 

40 



100 



165 



200 



100 



145 



0-1 



190 



- 2-4 
± 0-5 



{ 



170 
140 



4,000 



10,000 



100? 



200 \ 



2,000 



130 



1,400 



140 



100,000 



} 



340 



1270 



120 
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Eddy-stresses and of Eddy-diffusivity. 

the level of observation, Oy horizontally to the left and Oh vertically upwards. All 
are positive when the air above drags the air below in senses of x and y increasing. 



yy 

dynes per 
square 
centi- 
metre. 



hh 

dynes per 
square 
centi- 
metre. 



2*9 
0*6 



0-004 
0-006 



0*2 



-0'04 



-0-6 
±0-12 



0'7 



-0-34 



dynes per 
square 
centi- 
metre. 



yh 
dynes per 
square 
centi- 
metre. 



[ + 0-45] 
±0*39 



-0-48 
±0-20 



hx 
dynes per d v x 
square 



dh 



centi- 
metre. 



sec* 



0-1 



[-0-34] 
+ 0-18 



>0*5 



}0-3{ 



<0-06 



Conductivity 

'= pK 

grm. cm."" 1 sec." 1 , 



0-006 



0-07 



0-03 



0-3 

0-9 



12 



0-25 



6-8 



5-9 




6 




0-16 



Diffusivity. 
K 

cm. 2 sec."" 1 . 



240 

750 



10,000 



200 



5,500 



4,800 



5,000 



130 



Turbulivity. 

£ = g 2 p 2 K 

grm. 2 cm. 2 



sec. 



-5 



5 


7 


60 


80 


24 


34 



340 
1,050 



14,000 



300 



8,000 



7,000 



6,000 



200 
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Table IY. — Observed Values of 



Entropy gradient 
vertically upwards. 


Surface. 


! 

Methods. 


Notes. 


h 

cm. 

above 

surface. 


i 
cm. 

• 

sec. 


Ji/Ju 

dynes per 

square 
\ centi- 
! metre. 


1919, Oct., 21d. 12h., 
Benson 


Obstructions up 
wind -^ radian 


Paraffin vapour . . . 




160 


250 


— 


1919, Oct., 29d. 15h., 
Benson; raining 


Obstructions up 
wind • 045 
radian 


Smoke ..... 




190 
900 


330 


[ 


Small? 


i 

1 Moor, trees . . 


Burning rubbish . . 


VII. 


300 


': 


Small? 


Fields, trees . . 


■ 
Factory chimney . . 


VIII. 


3,000 


600 

1000 
200 


Negative ? . . . . 


Moor, low hills . 


Cloud ...... 


IX. 


150,000 
2,000 


— 


Small? 


Wooded hills . . 


Large fire ... . 


X. 


25,000 
37,000 


1600 






Shell-puffs .... 


XI. 


300,000 


500 


-80 


Positive ? . . . . 


Spurn head . . 


Dines anemometer 


XII. 


— 


1200 


Positive ? . . . . 


Open sea . . . 


Steamer's smoke . . 


XIII. 


1,500 


1 


Negative ? . . . . 


Open sea . . . 


Steamer's smoke . . 


XIV. 


5,000 


200 




- * . . 


Open sea . . . 


Alto-stratus . . 


XV. 


500,000 


— 





Notes to Table of Eddy-Stresses and Diffusivities. 



I. 



II. 



III. 



1917, June, 16d. 4h. 8m. L.A.T., hilltop near Ancemont, France. Standing hay composed of a 

species of Festuca (identified by my friend, Mr. Sam Pim). It grew fairly densely to 30 cm. from 

the ground and tall seed stems rose to 70 cm. 
1917, July, 16d. 19jh. L.A.T., Maffrecourt, France. Green oats 70 cm. high. No trees near. 

About sunset. Overcast with stratus. Observers : David Long and L. F. Bichakdson. 
1917, June, 29d. 19h. 30m. L.A.T., Viel Dampierre, France. A field of corn 60 cm. high. 

Clouds (stratus) motionless. Observers : F. H. Weatherall, G. Hutchinson, L.F.R. 
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Eddy-stresses and of Eddy-diffusivity (continued). 



i 

yy '1 

dynes perj 
square 
centi- 
metre. 


M 
dynes per 
square 
centi- 
metre. 


xy 
dynes per 
square 
centi- 
metre. 


yh. 
dynes per 
square 
centi- 
metre. 


hx 
dynes per 
square 
centi- 
metre. 


dh 
sec." 1 . 


Conductivity. 

- P K 

grm. cm." 1 sec. -1 . 


Diffusivity. 

K 
cm. 2 sec." 1 . 


Turbulivity. 

t = <?VK 

grm. 2 cm. 2 

sec."" 5 . 


1 
■~*~" i 

i 






. — . 


. — 


— 


2-5 


2,000 


3,000 




— . 


— . 


9-3 


7,500 


11,000 


ti | . ml 





— — 


, — 


, — 


6 


5,000 


7,000 


— - 


— 




12 


10,000 


14,000 


— 


-5 
















— 




— 




— 




120 


10 6 


l-4xl0 6 


— 


— 


— 


— 


. — 


— 


<0-9xl0 3 


<10 6 


<0-8xl0 6 


-110 






■ — 












— — 








— . 


0-12 


100 


140 








I 


— — 


12 


10,000 


14,000 


— 


-1 N 















IV. 



V. 



VI. 
VII. 



1917, June, 25d. 20h. 5m. L.A.T., Joinville, France. Moor with herbage dense to 10 cm. and 

rising thinly to 50 cm. 
1917, October, 4d. 8|h. G.M.T., Massiges, France. Flat moor with grass to 10 cm. and stems rising 

to 30 cm. Trees up-wind subtending an angle of 10 degrees. Overcast with strato nimbus, of 

which velocity/height = 0*025 sec." 1 . Temperature 0°*15 C. Observers : Olaf Stapledon and 

L. F. Eichardson. Eddies partly due to observer. 
See photographs and description in this paper. 
1917, July, 18d. 1\h. L.A.T., Maffrecourt, France. Moor, with small trees, 5m. high and houses. 

Overcast, 
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VIII. 1917, August, 2d. 17b. L.A.T., South England. Overcast. 
IX. 1917, July, 17d. 19h. 50m. L.A.T., East Champagne, France. Nine observations of a cloud at 

intervals of 1 minute by an Abney level. Angular elevation about 9 degrees. 
X. 1917, August, 8d. 20h. L.A.T., Argonne, France. Large fire of petrol and wood. Smoke observed 
from distance of 10 km. with a sextant. If smoke were not hot, height of its upper edge above 
ground would have given K = 4*6 x 10 5 cm; 2 sec. -1 , an over estimate. Irregularities in upper 
edge gave K = 10 4 , probably an under estimate. Mean K of order of 10 5 . Overcast. 
XI 1918, April, 12d. 14h. 5m. L.A.T., France. Two anti-aircraft shell puffs at a mean elevation of 
21 degrees above the horizon and 4 degrees apart in a vertical plane, were brought into 
coincidence in the field of view of a sextant. In 120 seconds their separation of 4 degrees did 
not vary visibly, certainly not by 4 minutes of arc. Their apparent motion was horizontal at 
tbVo' ra dian per sec. Height assumed 3200 metres — a likely value. 
XII. 1914, November, 16d. 3h. to 9h. G.M.T. Taken from the re-production of the Dines anemogram 
on p. 81 of the ' Observer's Handbook/ Meteorological Office, London, 1917 edition. 

XIII. 1917, July, 26d. 21h. L.A.T., the English Channel. 

XIV. 1917, July, 26d. 16h. L.A.T., the English Channel, off Havre. 

XV. 1917, July, 26d. 15h. L.A.T., the English Channel, off Havre. Observations of cloud at intervals 
of 1 minute with pocket sextant. Angular elevation about 4 degrees. 

The photograph (2) also shows that the smoke spreads more rapidly upwards than 

downwards, indicating that the stress hh and the turbulivity £ both increase with 
height. 

Fig. 3 shows another case of: low eddy-conductivity occurring at sunset. The 
smoke here is from burning hydrogen phosphides ; it is warm and rises slightly. To 
the eye the smoke appeared as a narrow wavy ribbon moving with a mean velocity of 
1'2 metres per second. The broader smooth band shown in the photograph is due to 
the exposure of 75 seconds, made long in order to get an average effect. The source 
of smoke is a bottle 3*4 metres above ground and just within the picture. The 
bamboos are 5 metres apart. The air density was 0*00126 c.g.s. The eddy- 
diffusivity works out to about 130 c.g.s. units, the eddy-conductivity to 0'16 c.g.s., 
the turbulivity to 200 c.g.s. The sky was cloudless. Obstructions to windward rose 
above the horizon to an angle of only §V radian. The photograph was taken in 
latitude 51° 37' N., longitude 4m. 24s. west, at 1919, Sept., 29d. 17h. 58m. G.M.T. 

Above is a table of observations. It is noticeable that when two of the 

direct stresses xx 9 yy> hh have been measured at the same time and place, they 
have been found to be not very unequal. G. I. Taylor has published some 
observations which show the same thing. It is as though there were a kind of 
equipartition of energy between the three components of the eddying motion. A 
very marked increase in both direct stress and diffusivity takes place either with 
velocity or with height. A rapid increase of viscosity with height in the first 
200 metres has also been deduced by W. Schmidt from wind observations made by 
Hellmann over a piece of flat land. (' Sitz ? Akad. Wiss./ Wien, 1917, Heft 6, 

p. i7)- 



Richardson. 



Phil, Trans., A, vol. 221, Plate 1. 
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The six observations of steamers' smoke are put forward only as upper limits to the 
turbulence appropriate to the bare sea, for the steamer itself probably makes 
a considerable eddy. 

It would be desirable to classify the observed eddy-conductivity as a function of 
four independent variables ; namely, the height, the vertical gradient of entropy, the 
vertical gradient of velocity and the character of the surface. Vertical gradient of 
velocity is suggested as an independent because it measures the only rate-of-mean 
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strain which attains a noticeable value in the free atmosphere, and because Osborne 
Reynolds* has shown that the energy of the eddy motion comes from the work 
done by the eddy-stresses upon the corresponding rates of mean strain. The observa- 
tions here presented are much too scanty for such a classification, but to render the 
relation to height visible, the effect of velocity has been removed, in one sense, by 
dividing each value of the eddy-conductivity by the velocity at that level. The 

* Lamb, ' Hydrodynamics/ IV. edition, § 369, equation (21). 
VOL. CCXXI.— A. E 



26 MR. LEWIS F. RICHARDSON ON 

justification for this procedure is that Taylor has given reasons^ for supposing that 
the viscosity, and therefore also the conductivity, is proportional to the velocity. 
For comparison with the present observations, the diagram shows Taylor's mean 
value of the diffusivity at the Eiffel Tower, f and also some general meansj deduced 
from precipitation by the writer. 

In order to compare them it has been necessary to assume some corresponding 
velocities; for which purpose I have taken 540 cm. sec." 1 at the mean height of the 
Eiffel Tower, 700 cm. sec." 1 as a world-mean at 500 metres and 1000 cm. sec."" 1 for 
the same at 8500 metres. These are based on information given in Hann's 
' Meteorology.' The conversion formulae between eddy -conductivity, diffusivity and £ 
have been given in Section I. In order to compress into a diagram the large ranges 
of height and conductivity, logarithms have been plotted. A smooth curve is drawn 
through the clustered observations over land. It shows a maximum between the 
heights of 100 and 1000 metres, and a marked falling off above and below. Not 
only c/v but also c the conductivity has a maximum here. Taylor's first 
observations related to heights near this maximum and so he naturally came to the 
conclusion that there was no marked variation with height. 

IX. Cumulus Eddies in Calm Weather. 

The familiar sequence, which can be observed in many places, is here illustrated by 
the mean of some selected days in latitude 49° in France, on the bare grass moors to 
the west of the forest of Argonne, in the month of May. The sun rose at 4h. 20m. 
local apparent time, but could not be seen for mist. By 6h. the disk of the sun 
became visible. At 7|-h. the mist was rising in large pieces, leaving a brilliant blue 
sky. At 9h. the first cumuli appeared over the forest. About half-an-hour later 
they appeared over the grass land also. By noon the cumuli covered -^ of the sky. 
By 16h. the cumuli had begun to spread out horizontally, and by 19h. they had 
vanished, leaving the sky clear again. 

Now here we have a collection of eddies in which the rising parts, represented by 
the cumuli, visibly move to a level where they remain by mixing with their 
surroundings. So we should be able to calculate the diffusivity K by the direct 
application of the formula given by G. I. Taylor (' Phil. Trans./ A, vol. 215, p. 3) 






aJj v K( h - hf ) dxd y> • i 1 ) 



where h—h! is the height through which the air has moved before mixing, v K is its- 
vertical velocity, and A is a large horizontal area. Only, as Taylor's formula assumes 

* G. I. Taylor, 'Boy. Soc. Proe./ A, vol. 92, pp. 196-199. 
t G. I. Taylor, 'Roy. Soc. Proc.,' A, vol. 94 (1917), p. 141. 
J L. F. Richardson, ' Roy. Soc. Proc./ A, vol. 96 (1919), p. 18. 
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that p and K are independent of height, it may be as well to remove these restrictions. 
It is then found that £ as defined by equation (l), of Part L, is given by 

£ = 2i m ^ (p f ~-p) dx dy, • (2) 

where p f is the pressure at the initial level, m H the vertical momentum per volume. 

Now let us insert numerical values. The surface air was seen to begin to move at 
or a little before 7^h. The cumuli appeared to cease rising before 16h. and the height 
of their tops is known to average about 2 km. (vide Hann's ' Meteorology/ Illrd edn., 
p. 280). Now if we suppose, as seems reasonable, that the top of the cumulus is formed 
from the damp air which was initially close to the ground, then the displacement, 
measured by pressure, is about 2 decibars, so that (p'—p) = 2 x 10 5 dyne cm." 2 . The 
vertical velocity is 2 km. in 8*5 hours, that is 6*5 cm. sec." 1 . So the momentum per 
volume = m — 7 x 10" 3 grm. cm." 2 sec. -1 in the rising current. The rising current 
covered 0*4 of the sky, so that averaging over the area A, as is done in (2), is 
equivalent to taking 0*4 of m K (p'—p) for the rising current. But the invisible 
descending currents contribute an equal amount to the integral. So 

£ = 0xO , 8x2xlO 5 x7xlO- 3 
= 11 x 10 5 grm. 2 cm." 2 sec." 5 . 

This figure is about ten times greater than measures of £ at a height of a few 
hundred meters, deduced by various authors. If the air which forms the top of the 
cumulus had really started from a height of 1 km. instead of from the ground, as we 
have supposed, then the numerical value of £ would have to be divided by four. 

Keasons have already been given (Part VII.) for supposing that £ derived in this 
way from cumulus clouds is a measure of frictional effects, but not of the diffusion of 
entropy, because the linear term in (Part V., 3) does not vanish on taking the mean, 
owing to the fact that the eddies are produced by variations of entropy. To put it in 
another way : In G. I. Taylor's deduction of formula (l) the vertical gradient of the 
diffusing quantity is treated as not correlated with the vertical velocity. When we 
are dealing with cumulus clouds that assumption is probably justified if the diffusing 
quantity is horizontal velocity, but not if it is potential temperature. 

To find £ in the sense of diffusivity for potential temperature we should have to 
employ formula 32 of Part VII., namely 



K 



/39V f 



I — 1 ^H 



d 2 e/dtt 



For insertion in this we require lapse rates in cumulus clouds and in the clear air 
between them. Such have recently been obtained by airmen. 
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X. Summary. 

Part I. deals with notation. Measures of turbulence may advantageously be 
expressed in the form (p in 

3 X _ 9 (* 9 X 



dt dp \ dp 

where p is the pressure (here used as a measure of height), t the time, and x may 
be horizontal velocity in a fixed azimuth, or potential temperature, or water per 
mass of atmosphere. It is suggested that £ might be called the "turbulivity." Its 
dimensions are grm. 2 cm." 2 see." 5 . Better still is the conductivity c = £p~ l g~ 2 - 

In Part II. the eddy-shearing stress on the ground is deduced from pilot balloon 
observations. Values on land in any self-consistent dynamical units are found to 
range from 0'0007 to 0'007 times the value of m 2 /p, where m is the mean momentum 
per volume up to a height of 2 km. and p is the density. Compare G. I. Taylor, 
6 Roy. Soc. Proc.,' A, vol. 92. 

In Part III. evidence is given to show that the eddy-viscosity across the wind at 
Lindenberg increases with height, and, except near the ground, is much greater than 
the eddy-viscosity along the wind. Here £ ranges from 10 4 to 5 x 10 5 . 

In Part IV. the spreading of a lamina of smoke is considered. Values of £ ranging 
from 7 to 140,000 are found. £ increases both with height and with velocity. 

In Part V. the derivation of £ from smoke observations is examined more 
thoroughly. 

Part VI. deals with Osborne Reynolds' eddy-stresses. For one occasion an 
attempt was made to measure simultaneously all six components of stress by 
observing .the motion of thistledown. The three direct stresses are easily measured. 
Not so the shearing stresses however, one was found to be 2*4 times its probable 
error. 

Part VII. summarizes the theory of scattering of particles. 

Part VIII. contains numerical values derived from scattering. 

In Part IX. the turbulivity £ is estimated from the rising of cumuli in calm 
weather and found to be 10 6 , applicable only in the sense of friction. Thus the 
whole range of £ observed in the free atmosphere was from 7 to a million in 
contrast with 0*2 in perfectly still air in a laboratory. The eddy-stresses observed 
have ranged in absolute value from 0*004 to 110 dynes cm. 
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